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^ ; ABSTRACT 

^ '. We use SDSS data to investigate the scaling relations of 127 NoSOCS and 56 CIRS 

Qh' galaxy clusters at low redshift (z ^ 0.10). We show that richness and both optical and 

q X-ray luminosities are reliable mass proxies. The scatter in mass at fixed observable is 

;_i ' ~ 40%, depending on the aperture, sample and observable considered. For example, for 

' the massive CIRS systems tr/ n M500|Af500 = 0.33 ± 0.05 and <Ji n M500\Lx = 0.48 ± 0.06. 

C3 . For the full sample 07„M5OO|iV5OO = 0.43 ± 0.03 and cri n M500\Lx = 0.56 ± 0.06. The 

scaling relations based only on the richer systems (CIRS) are slightly flatter than those 
I \ based on the full sample, but the discrepancies are within l-a. We estimate substruc- 

■ ture using two and three dimensional optical data, verifying that substructure has no 
(*T) ' significant effect on the cluster scaling relations (intercepts and slopes), independent of 
I/"") \ which substructure test we use. For a subset of twenty-one clusters, we estimate masses 
r***- ■ from the M-T^ relation using temperature measures from BAX. The scaling relations 
C*~) ' derived from the optical and X-ray masses are indeed very similar, indicating that our 

, method consistently estimates the cluster mass and yields equivalent results regardless 

■ of the wavelength from which we measure mass. For massive systems, we represent the 
| mass-richness relation by a function with the form ln(M 2 oo) = A + B x ln(N 2 oo/60), 

with M200 being expressed in units of 10 14 M Q . Using the virial mass, for CIRS clus- 
ters, we find A = (1.39 ± 0.07) and B = (1.00 ± 0.11). For the same sample, but 
• i-H [ using the masses obtained by the caustic method, we get A = (0.64 ± 0.14) and B 

■ = (1.35 ± 0.34). If we consider the mass as estimated from Tx (for the subset of 21 

clusters with T x available) we derive A = (0.90 ± 0.10) and B = (0.92 ± 0.10). The 
relations based on the virial mass have a scatter of cri n M20Q\N200 = 0.37 ± 0.05, while 
f7 inM200|A f 200 = 0.77 ± 0.22 for the caustic mass and cr;„ M2 oo|Af20o — 0.34 ± 0.08 for 
the temperature based mass. 
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1 INTRODUCTION 

Two aspects of galaxy clusters make their study uniquely 
important for cosmology. The path from density pertur- 
bations to virialized objects involves the physical mecha- 
nisms that ultimately determine the cluster properties as 
we currently observe them. The growth of these systems 
depends on cosmological parameters like f2 m , erg, etc. 
The sensitivity to these parame ters depends on the redshift 
jCarroll. Press, fc Turned 1 19921 ). Therefore, the systematic 
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study of cluster properties, based on unbiased samples at 
different redshifts, can provide substantial input for models 
trying to explain not only how clusters formed and evolved 
from the astrop hysical viewpoint but also how spacetime is 
structured (e.g. iHenrv fc Arnaudl fl99ll ; iRosat et all 120021 ). 
Ultimately, we aim to disentangle the astrophysical aspects 
from the cosmological ones. 

Clusters have no sharp edges. Observations show that 
their central parts reached virial equilibrium earlier while 
the outskirts are still accreting matter from their surround- 
ings, making the measurements of cluster properties a dif- 
ficult issue. Clusters are constituted of both dark (non- 
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baryonic) and luminous (baryonic) matter. This mixture 
renders observations within a specific wavelength range very 
limited in utility and determining the physical state of 
each component too complex. Over the years, the sizes of 
cluster samples and the quality of the data hampered our 
progress on all of the fronts mentioned above. Only re- 
cently have large samples with high-quality data been ob- 
tain ed, allowing a systematic investiga tion of these effects 
(e.g. IVikhlinin et alll200Sl ; lHenrvll2004r i . 

The abundance of clusters of gal axies can be us ed as 
a probe of linear density fluctuations (|Eke e t al. 1996j) and 
since their growth is directly determined from the dark en- 
ergy properties, measuring the number of clusters of a given 
mass as a function of redshift can constrain the equation 
of state (w p arameter) of this majority compo nent of the 
U niverse (e.g. [ Frieman. Turner, fc Hutererll200ll ). As shown 
bv lHuterer fc Turner! l|200ll) . the sensitivity of cluster abun- 
dance to w is maximum around redshift 1.0, so that a cluster 
sample over a large area on the sky and deep (z ~ 1.0) is 
critically needed to provide reliable results. Direct measures 
of halo mass are not possible, so that we have to rely on 
mass proxies free of systematics. The most accurate mass 
tracers are velocity dispersion and X-ray temperature. Un- 
fortunately, for large samples, at any redshift, directly mea- 
suring these cluster properties is not feasible. Hence, we have 
to rely on other observables (such as richness or X-ray lu- 
minosity) that are easier to derive for all clusters in a given 
survey. 

Several studies have investigated mass- observable (MO) 
relatio ns and their evolution with redshift in the optical and 
X-ray (IVikhlinin et al.ll2002l ; lLevine et aHl2002l ; lEttori et al ' 



2004 ; IStanek et alj 12006); iBecker et all 120071; iRykoff et al 
200a ). As advocated by [Lima fc Hul l|2004 120051 ) the 



ance of the counts themselves (clustering of clusters) can 
be used to calibrate the MO relation. Dark energy depends 
only on the redshift, but the cluster properties vary with 
mass and redshift. Using cosmological simulations one can 
use counts and their variance to not only normalize the MO 
relation but also learn which cosmology fits the data bet- 
ter. Therefore, understanding the scaling relations of galaxy 
clusters and their dependence on the dynamical state of the 
system is crucial. 

This paper focuses on four specific issues: (i) examin- 
ing the performance of different cluster properties, like rich- 
ness (Ng a i s ), optical (L opt ) and X-ray (Lx) luminositief], as 
proxies for cluster mass. The relation between one of these 
proxies and an independent mass estimate defines a scaling 
relation; (ii) measuring how much substructure there is in 
clusters (in 2D and 3D) and establishing its effect on the final 
mass calibration; (iii) investigating how the scaling relations 
change when mass is derived from the analysis of the veloc- 
ity distribution or temperature. Optical masses used in this 
work are eithe r from a virial analysis ( ours) or from the caus- 
tic technique l|Rines fc Diaferioll2006l . hereafter RD06); and 
(iv) comparing the optical and X-ray properties of galaxy 



1 The error in the X-ray luminosity presented in the first paper of 
this series was underestimated. That is fixed in the current paper. 
All the X-ray measurements available in Table 5 of the first paper 
are replaced by the values listed in Table Al of this paper. 



clusters, as well as obtaining a direct relation between R200 
and richness. 

This paper is organized as follows: §2 briefly describes 
the samples used here, while §3 presents the mass calibra- 
tion based on N ga i s , L opt and Lx- In this section we also 
investigate the effects of substructure on the scaling rela- 
tions and compare the results based on different mass esti- 
mators. These are derived with optical data, from the virial 
or caustic analysis, or from the gas temperature. The mass- 
to-light ratio is also presented at the end of this section. In 
§4 we discuss how our results compare to others in the lit- 
erature. Correlations between optical and X-ray properties 
are presented in §5, while in §6 we establish the connection 
between richness and physical radius. The main conclusions 
are drawn in §7. The cosmology assumed in this work is 
fi m =0.3, fl x =0.7, and H = 100 h km s" 1 Mpc" 1 , with h 
set to 0.7. 



2 DATA AND METHODS 

The main data set used in this work is the supplemen- 
tal version of the Nor t hern Sky Optical Cluster Survey 
(NoSOCS, lLopesI 120031 : lLopes et alll2004h . which has its 
origin on the digitized version of the Sec ond Palomar Ob- 
serva tory Sky Survey (PO SS-II : DPOSS. IpTorgovski et ail 
l2003h . In iGal et all (|2004l ) and lOdewahn et all (|2004l ). we 
describe photometric calibration and object classification for 
DPOSS, respective ly. The supplemental version of NoSOCS 
|Lopes et all 120041 ) goes deeper (z ~ 0.5), but covers a 
smaller reg ion than the main NoSOCS catalog l|Gal et all 
2003, 2009), and contains 9,956 cluster candidates over ~ 
2,700 square degrees. The smaller area is due to the use 
of the best DPOSS plates, selected according to seeing and 
limiting magnitude (r = 21.0). 

We examine a sample of low redshift galaxy clusters 
(z ^ 0.10) from the NoSOCS supplemental catalog. As this 
survey comprises mainly poor systems in this redshift range, 
we complemented them with more massive systems from 
the Cluster Infall Regions in SDSS (CIRS) sample (RD06). 
CIRS is a collection of z ^ 0.10 X-ray selected clusters over- 
lapping the SDSS D R4 footprint. As d escribed in the first 
paper of this series l|Lopes et all 12009) we extracted SDSS 
data for 127 NoSOCS clusters and 56 CIRS systems at low- 



We used photometric and spectroscopic data from 
the fifth release (DR5) of the Sloan Digital Sky Survey 
<| York et al.ll2000h . The exception is for the CIRS systems, 
which were incorporated in a late stage of this work, thus 
having data from the SDSS DR6. All the magnitudes re- 
trieved from SDSS are de-reddened model magnitudes. De- 
tails r egarding SDSS da ta extraction are provided in pa- 
per I (|Lopes et all 120091 ). There, we used SDSS photomet- 
ric dat a to estimat e more accurate photometric redshifts 
(Zphoto, lLopesll2007M . richnesses and optical luminosities for 
the full NoSOCS supplemental catalog. We found 7,414 sys- 
tems well sampled in SDSS DR5. Appr oximately 10% (754) 
have z p hoto ^ 0.133 (|Lopes et alll2009T ). 

The above redshift limit comes from the choice of only 
using clusters at z ^ 0.10 for the current work. In this red- 
shift range t he SDSS spectrosc opic survey is complete. As 
discussed in lLopes et all l|2009l ). at higher redshifts galax- 
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ies fainter than M* + 1 are missed, bia sing the dynamica l 
analysis (see discussion in section 4.3 of lLopes et al.ll2009l ). 
Out of the 754 NoSOCS supplemental clusters with z p hoto ^ 
0.133, we were able to determine the spectroscopic redshift 
for 179 systems, requiring at least 3 galaxies within 0.50 h~ 
Mpc in the SDSS spectroscopic footprint. Note that we only 
selected systems for which \z p hoto — z spec \ < 0.03(1 + z spec ). 
We eli minated interloper s using the "shifting gapper" tech- 
nique (|Fadda et al.lll996r i. applied to all galaxies with spec- 
tra available within 2.50 h _1 Mpc. From the 179 clusters 
with z apec ^ 0.10 we retained 127 systems with at least 
10 member galaxies selected by the above procedure. We 
applied the sa me procedure to th e 56 CIRS clusters. Fig- 
ures 4 and 5 of lLopes et ail (|2009h show the velocity-radius 
distributions of the 127 NoSOCS clusters and the 56 CIRS 
systems. 

These clusters were then subjec ted to a virial analy- 



sis analogous t o the one d escrib ed in Girardi et al.l (I X99; 



IPopesso et all (|2005l . 120071 ) and iBiviano et ail |2006h . This 
procedure yields estimates of ap, -R500, -R200, A/500 and M200 
for the set of 183 low redshift clusters considered in this 
work. X-ray luminosity is also estimated for these systems 
using ROSAT All Sky Survey (RASS) data. A detailed de- 
scription of both the vir ial analysis and de rivation of X-ray 
luminosity is provided in lLopes et all l|2009h . which also lists 
the properties of these systems in Tables 1, 2, 3, 4, and 5. 

The NoSOCS clusters have velocity dispersion esti- 
mates of 100 < ap < 700 km/s. The CIRS systems have 
200 < a P < 900 km/s (with only 23% of objects with 
ap < 400 km/s). For most of the studies below we ignore 
three of the CIRS objects which have biased values of ap 
and mass due to projection effects and substructure (see pa- 
per I). We only show these objects in Figures 9, 10, and 11, 
where we investigate the impact of substructure on the scal- 
ing relations. These 3 clusters are Abell 1035B, Abell 1291A 
and Abell 1291B. 



3 MASS CALIBRATION 

Scaling relations involving a simple cluster observable, like 
richness or luminosity and a fundamental property such as 
mass, can provide important clues on how large scale struc- 
ture forms, galaxy formation proceeds, and how the intra- 
cluster gas reaches its current state. The establishment of 
these relations for nearby clusters is also vital for future 
studies of the cluster population in the distant universe. Par- 
ticularly, th^y_jna^_liel_rj_£c^trahi the dark energy equation 
of state jMaiumdar fc Mohrlliooi 'l. 

In this section, we examine the correlation of N ga is , L opt 
and Lx with cluster properties, such as velocity dispersion 
(ap) and mass (M200 and Msoo). We use the spectroscopic 
redshifts of the clusters (instead of z p hoto) when comput- 
ing richnesses and luminosities. Three different apertures are 
used to derive these quantities: 0.50 h _1 Mpc, R500 and R2oo- 
We consider the re-centered (luminosity-weighted) coordi- 
nates for the 127 NoSOCS clusters and the original position 
(X-ray centroid) listed in RD06 for the 56 CIRS clusters. 
We present the results for the NoSOCS and CIRS samples 
independently as well as for the combined sample. 

In the next three subsections we discuss the mass cal- 
ibration based independently on richness, optical luminos- 



ity, and X-ray luminosity. For the richness relations we also 
investigate the impact of using less rigorous criteria for in- 
terloper removal when computing mass, as well as assuming 
different centroids and considering all NoSOCS clusters at 
2 ^ 0.25. We then examine the impact of cluster substruc- 
ture on the scaling relations and compare the results derived 
from mass estimates obtained from different wavelengths. 
Finally, we show the dependence of the mass-to-luminosity 
ratio with the cluster scale (defined by its mass). 



3.1 Calibration with richness 

Figure 1 shows the comparison between mass (M500 and 
M200) with richness, estimated within R500 and R2oo- The 
lower panels show the residuals, LOG(M„i, s /M/ji). M b s is 
the observed value of the mass within the given radius, 
while M/a is the linear regression solution. The results 
are only shown for the 127 NoSOCS objects. The solid 
line in each upper panel sho ws the orthogonal regression 
fit (|Akritas fc Bershadvlll99"6l ). All the scaling relations ob- 
tained in this work are of the form 

ln(Y) = A + Bxln(|) , (1) 

The X and Y parameters are listed in the tables defining the 
relations (see below). The pivot point (C) depends on the 
sample being used. It is taken to be approximately the me- 
dian value of the X parameter within R500 for a given sam- 
ple. For the NoSOCS and full samples (NoSOCS + CIRS) 
we consider the same pivot point, while for the more mas- 
sive systems (CIRS) we have another pivot point. Those are 
listed in the tables and can be seen in the figures. The linear 
fit is obtained by a two iteration 3-sigma clipping, so that 
after a first run we eliminate outliers. Those are not used in 
the final fit and are indicated (if exist) as open symbols in 
the figures. 

The results for the scaling relations involving ap, M500, 
M200 and N ga ! S are summarized in Table 1. The columns give 
(1) the parameters involved in the relation (the abscissa and 
ordinate of the corresponding plots); (2) the cluster sample 
used; (3) and (4) the intercept (A) and slope (B); (5) the 
scatter (in natural log space) in the Y parameter at fixed X; 
(6) the total number of clusters used in the preliminary fit; 
and (7) the final number after a 3-a clip. Note that we eval- 
uated the error in the scatter, performing a boostrap proce- 
dure (500 events) with replacement. So the scatter is shown 
with the associated error in all tables (exception to Table 8). 
Rows 1-6 provide the results using only the NoSOCS clus- 
ters, as shown in Figure 1. Rows 7-12 are for clusters from 
RD06 only, while the combination of NoSOCS and RD06 is 
listed in rows 13-18. 

If one has to work within a fixed metric radius, 0.50 
h _1 Mpc yields the smallest error in the richness or L opt 
measures and minim izes the scatter in the optical versus X - 
ray scaling relations (|Popesso et al ]|2004l ; lLopes" et al 1120061 ). 
However, when compared to apertures that scale with mass, 
this radius (0.50 h _1 Mpc) is too large for low mass sys- 
tems and too small for massive ones, and may lead to tilted 
(steeper) scaling relations (see Table 1). When we use an 
aperture that scales with mass the relations become flat- 
ter. For instance, for the NoSOCS sample, the slope of the 
M2oo-N 9a ; s relation is 1.45 when using the fixed metric ra- 
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Figure 1. The connection between mass and richness, computed with two apertures (R500 and R,20o)- In the lower panels the residuals 
LOG(M D b s /Mji t ) are shown. M D £, S is the observed value of the mass, while Mj it is the linear regression result. 



dius for computing richness (0.50 h 1 Mpc), but only 1.09 
when using R.200- 

From the lower panels of Figure 1 we see that some 
of the poorest systems (N^° < 30) have larger devia- 
tions from the linear solution. Both richness and mass are 
harder to determine due to the low contrast of some of these 
clusters. Incompleteness in spectroscopic sampling for these 
poor systems may also be a problem (see discussion in §3.2). 

The comparison of the relations for the 127 NoSOCS 
clusters, the 53 CIRS systems and the combined sample, 
shows that the scatter of the relations obtained with the 
CIRS and extended samples is smaller than for the NoSOCS 
clusters alone. This is because CIRS is restricted to more 
massive systems and the enlarged sample has better statis- 
tics, spanning a wider range of all parameters. For the M200- 
~^gais relation the orthogonal scatter in mass at fixed rich- 
ness is (0.60 ± 0.07) for NoSOCS, (0.37 ± 0.05) for the CIRS 



sample and (0.46 ± 0.04) for the combined data set. The 
slopes are nearly the same for the NoSOCS and extended 
samples, but slightly different for CIRS, although they are 
compatible within the errors. For the M2oo-N^ i ° relation 
the slope is 1.09 for the NoSOCS sample, 1.00 for CIRS and 
1.07 for the combined sample. 



3.1.1 Systematics in the scaling relations 

Table 2 lists the coefficients of the relations ap-N^^ and 
M200-N™" using only the NoSOCS clusters, with three dif- 
ferent methods. In the first two rows we consider the original 
cluster coordinates instead of the luminosity- weighted ones. 
Rows 3-4 list the results when we use a more "relaxed" cr i- 
teria for interloper removal (see § 4.1 of lLopes et al.ll2009h . 
and in rows 5-6 we consider all NoSOCS clusters at z 
0.25, instead of z ^ 0.10 (limit adopted for this work). The 
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latter is useful to see how much the scaling relations are af- 
fected for using clusters at redshifts where the spectroscopic 
sampling is ve r y inco mplete (see § 4.3 of paper I) . Note that 
iPopesso et alj (|2005t ) (hereafter POP05) considered clusters 
at z ^ 0.25 for deriving scaling relations based on SDSS 
data. 

When comparing the results obtained with the original 
coordinates (first two rows of Table 2) to the ones with the 
luminosity- weighted centroids (Table 1), we find that the 
intercept, slope and scatter are consistent within the errors 
(1-cr). Thus, we conclude that the centroid is not a critical 
issue when deriving the scaling relations. When the criteria 
to select interlopers are relaxed (results in rows 3-4 of Table 

2) the agreement is not as good for the intercept and slope, 
but it is still within 1.5-a. Relaxing the criteria for rejecting 
interlopers in the "shifting-gapper" technique has a minor 
effect in the relations. 

The third case in Table 2 uses all NoSOCS clusters at 
z ^ 0.25 (rows 5-6). A comparison to Table 1 reveals that 
the intercept and slope are now different, in some cases con- 
sistent only within 2-cr. The scatter is, however, similar. This 
result shows the relevance of using complete spectroscopic 
samples for studying the scaling relations. The use of clusters 
at redshifts where the spectroscopic survey is not complete 
to M* + 1 results in biased velocity dispersion and mass 
estimates (see paper I), yielding flatter scaling relations. 

This is a critical issue for measuring the cluster mass 
function since the mass calibration might be severely biased. 
However, it is important to note that this bias in the calibra- 
tion may reflect the percentage of clusters with incomplete 
spectroscopic sampling. In this work 127 of the 219 NoSOCS 
clusters at z ^ 0.25 are below z ^ 0.10, and are therefore 
well sampled. Ninety-two of the 219 systems (42%) have 
poor sampling (z > 0.10), a considerable fraction of the to- 
tal used when deriving the scaling relations. Other authors 
included z > 0.10 clusters in their sample (POP05, for in- 
stance), but show results consistent to our unbiased sample 
(at z 0.10). Although they initially have many systems at 
2 > 0.10, they require that clusters have at least ten galax- 
ies with redshifts when estimating mass and deriving the 
scaling relations, which may preferentially exclude many of 
the higher redshift clusters. We believe that this is the case 
since their results are consistent with ours. 

3.2 Calibration with optical luminosity 

Figure 2 shows the comparison between mass (M500 and 
M200 ) and optical luminosity computed within R500 or R200 • 
Residuals are shown in the lower panels. We only show the 
results for the 127 NoSOCS objects. One interesting feature 
from the comparison of Figures 1 and 2 (see Tables 1 and 

3) is that the slopes obtained in the two cases are consis- 
tent. This indicates the optical luminosity and the number of 
gala xies used to compute L op t, are proportional to each other 
(see lPopesso et al1l2007h . with a constant of proportionality 
consistent with unity. Figure 3 shows op versus L op t, com- 
puted within R500 or R200 for only the CIRS objects. Figure 
4 is analogous to Figures 2 and 3, but shows the extended 
sample of 180 NoSOCS plus CIRS clusters. The inspection 
of the lower panels of Figures 2 and 4 reveals that larger 
deviations from the linear regression occur for L^ t 00 < 0.30 
10 12 L©, corresponding to the cut at N^j, — 30 seen in 



Figure 1. Hence, we stress the conclusion that richness, L opt 
and mass are harder to obtain due to the low contrast of 
such systems. Another problem affecting these clusters may 
be incompleteness in the spectroscopic sampling. 

To confirm these conclusions we separate the clusters 
shown in Figure 4 in systems with |LOG(M b s /Mfl t ) | ^ 2 
or |LOG(M ob s/M flt )| > 2. We find 135 clusters in the 
first case and only forty-five (25%) in the second. For 
L0G(M o b s /Mflt) I ^ 2 we find the following median val- 
ues: N 9ais = 33.0, N gais _ err /N 9ais = 0.24, N apec = 48.0 
and N spec 2oo = 21.0; where N ga ; s is computed within R200, 
N ga ; 3 -err/N 9a ; s is the percentage error on the richness mea- 
surement, N spec is the number of members (spectroscopi- 
cally selected) within the maximum aperture (normally 2.5 
h _1 Mpc; see paper I) and N spec 2oo is the number of mem- 
bers within R200. For |LOG(M obs /M fit )| > 2 we have N gais 
= 22.0, Ngcis-err/Ngais = 0.36, N spec = 30.0 and N spec2 oo 
= 15.0 Thus, the forty-five clusters with larger mass devi- 
ations are generally poor, low contrast (indicated by their 
large percentage error in richness) and have their mass esti- 
mates based on few galaxies. 

The results shown in Figures 2, 3 and 4 are summarized 
in Table 3. The meaning of all columns and rows is analo- 
gous to those in Table 1. The same conclusions drawn for 
Table 1 are valid now. The most important result is that the 
relations involving L op t have approximately the same scatter 
as those obtained with N ga i s . In other words, optical lumi- 
nosity performs as good as richness for mass calibration. For 
the full sample, the scatter of M200 at fixed N^ is (0.46 ± 
0.04), while it is (0.49 ± 0.04) at fixed L^ 00 . For the CIRS 
sample, the scatter of M200 at fixed N^ is (0.37 ± 0.05), 
being (0.36 ± 0.06) at fixed L^ 00 . 

3.3 Calibration with X-ray luminosity 

Next we compar e the X-ray luminosity estimated in 
lLopes et all (|2009l ) to ap and mass. As described in paper I, 
we computed the X-ray luminositj0, from RASS, using three 
different background estimates, termed "annulus" (from a 
ring surround the cluster), "boxes" (100 randomly selected 
background boxes) and "frame" (from the whole frame in 
which the cluster is located). Here we show the results ob- 
tained with the "annulus" and "frame" backgrounds. 

Figure 5 shows the relation between Lx ("annulus" 
background) and ap (top panels) and mass (lower pan- 
els) , estimated within R500 (left) and R200 (right). We lose 
seventy-four of the 180 clusters from the extended sample, as 
their X-ray luminosities are upper limits. These are mostly 
poor systems. When establishing the scaling relations in- 
volving Lx, upper limits were not used in the fits. From 
this plot (and Table 4) we see that the smallest scatter is 
generally found when using Rsoo- 

2 Due to a bug in the code used to estimate the X-ray luminosity 
and its associated error, the background was not properly deter- 
mined in paper I. That resulted in slightly different values for 
Lx- However, the errors of Lx were severely underestimated. We 
fixed the code and updated the values of Lx and its error (for 
all the three backgrounds) in the current paper. The information 
given in Table 5 of paper I is now replaced by the data available 
in Table Al, in the appendix of the present paper. 
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Figure 2. Analogous to Figure 1, but for mass and optical luminosity. 



We also studied the relations considering the CIRS sub- 
sample, with the comparison of M200 and L x 200 shown in 
Figure 6. It is particularly important to see whether the 
slope and mainly the scatter are reduced for this sample 
based on richer clusters. As found for N ga ( 3 and L op t, the 
slope also becomes shallower for the richer systems when 
considering Lx ■ Note that the slopes are consistent with the 
full sample (Figure 5 and Table 4) within 1-a. The slope of 
the M 2 oo-Lf 200 relation is (0.67 ± 0.06) for all clusters and 
(0.52 ± 0.08) for CIRS (see Table 4). The scatters in the re- 
lations based on the richer clusters are slightly smaller than 
those of the whole sample, but still within 1-a. The relations 
based on the "frame" background are similar to the ones ob- 
tained with the "annulus" , being the agreement within 1-a. 
Figures 7 and 8 show the results with the "frame" back- 
ground, for the full and CIRS samples, respectively. 

Tables 4 and 5 summarize the relations involving Lx- 
Table 4 gives the results for the "annulus" background, while 



Table 5 is for the "frame" background. As in the previous 
tables, the first column lists the two parameters being com- 
pared. The background type is listed in the second column. 
The sample used is in the third column, while the remaining 
columns are analogous to Table 1. In both tables we list the 
results for the full sample and for the CIRS systems. 

When comparing the performance of optical parame- 
ters (richness or L opt ) to Lx as a trace of the cluster mass, 
we find that the former perform a little better than the lat- 
ter. However, in the worst cases, the scatter is consistent 
within 1.5-a (see Tables 1, 3, 4 and 5). The main reason for 
these minor differences may lie in the fact that RASS data 
is shallower than the optical data from SDSS, so that Lx 
is determined less accurately than N ga j s and L op t. It is also 
important to stress that the optical properties may be used 
as reliable mass proxies, in the same way as Lx- This find- 
ing corroborates the results from POP05, who found that 
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Figure 3. Analogous to Figure 2, but using only the 53 clusters from RD06 and considering <rp and L op t. 



L op t gives slightly more accurate results than Lx (also us- 
ing SDSS and RASS data). 



3.4 Impact of substructure in the scaling relations 

In lLopes et al.l l|20091i we described the use of photomet- 
ric and spectroscopic data for NoSOCS and CIRS clusters 
to estimate the fraction of systems with substructure. Two 
specific tests have proven to be very sensitive to the presence 
of s uch disturbances in t he galaxy distribution inside a clus- 
ter JPinknev et alJll996h. T he first test is the DS, or A test 
l|Dressler fc Shectmanl 19881 ). which is a three dimensional 
test. The second is a two dime nsional test , called the sym- 
metry or P test, introduced by I West et al.l lll988h. Detailed 
descri ption s of both tests can be found in Pinknev et al.l 
l|l996t l and lLopes et al.l (|2009T i. For both tests, the signifi- 
cance level is determined with the aid of Monte Carlo sim- 
ulations. We set our significance threshold at 5%. 



The substructure tests are only applied to clusters with 
at least five galaxies within the aperture being considered. In 
paper I we showed that when using an aperture of R200 the 
rate of clusters showing significant signs of substructure is 
~ 21% for both the A and /? tests. The latter can be applied 
to galaxies with spectra or to the photometric data alone, 
while the former test requires redshifts. If we apply the J3 
test to the galaxies within a fixed metric aperture (1.5 h _1 
Mpc) the r ate of clusters with substructure rises to ~ 35%, 
in line with lLopes et al] (|2006h . who estimated substructure 
from two dimensional optical data for > 10,000 clusters. 

In this work we do not intend to investigate substruc- 
ture for each cluster individually. Our goal is to check the 
possible effect of substructure on the scaling relations con- 
necting M200 to Ngaij, and L opt . Four different cases are con- 
sidered. First we apply the A test to the 170 clusters with at 
least five galaxies with z spcc within R2oo- Second, we apply 
the /3 test to the same data set. Third, the /3 test is applied 
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Figure 4. Analogous to Figure 2, but for the combined sample of 127 NoSOCS clusters and 53 systems from RD06. 



to the photometric galaxy distribution within R 200 ■ We con- 
sider all galaxies with m* — 1 ^ m r ^m* + 1 fsee lLopes et al.l 
2006 for details) with no restriction regarding z ap e C - Fourth, 
we apply the f3 test to all galaxies within 1.5 h _1 Mpc of 
each cluster (the same photometric range is enforced). 

The results for the cases listed above are summarized 
in Table 6. In this table we always show first the results 
for the full sample and then the ones for clusters without 
substructure. The first six lines give the results considering 
the first two cases mentioned above. Lines 7 to 10 list the 
results for the third case, and the last four lines summarize 
the results for the fourth case. The first column lists the two 
parameters of the scaling relation (and the aperture adopted 
for computing N 9a ; s and L opt ). The second column shows 
the substructure test (if any) used to remove clusters with 
substructure. Note that the aperture used for estimating 
substructure is R200 in the first three cases and is 1.5 h~ 
Mpc for the fourth. The sample used is listed in column 3. 



The meaning of the remaining columns is the same as in 
Table 3. 

Figures 9, 10, and 11 show the results for the first, third 
and fourth cases discussed above. In each figure the relation 
between ap (top) and M200 (bottom) to L op t is displayed 
for all clusters (left) and for the substructure-free systems 
(right). From Table 6 and Figures 9 and 10 we see that 
the fitting parameters and scatter of the scaling relations, 
depend very little on the presence of substructure. The in- 
tercepts and slopes are always consistent within 1-a. It is 
important to note that the 2D and 3D tests indicate similar 
fractions of clusters with substructure, and the scaling rela- 
tions obtained for the substructure-free systems are nearly 
the same as those for the full sample, regardless of whether 
the test used to detect substructure is 2D or 3D. 

From Figure 11 we see that even considering the full 
photometric data within a fixed radius, instead of only clus- 
ter members within R200, the results are qualitatively the 
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Figure 5. The correlation between up, mass, and X-ray luminosity (estimated with the "annulus" background). The 76 clusters with 
upper limits (UL) for Lx are not included in the plot or used for the fitting. 



same. Due to the different apertures the scaling relations 
are different from the two previous figures. However, we still 
find the relations to be insensitive to the exclusion of clus- 
ters with substructure. The scatter of the relations (with or 
without substructure) shows a larger difference than for the 
other cases above, but it is still wi t hin 1 -a. These results 
contradict Figure 16 of lLopes et al.l ||2006T ). where a strong 
segregation in the Tx-N ffa ; s relation was caused by the pres- 
ence of substructure in some clusters. Several reasons may 
be responsible: 1) Tx was taken from the literature, so the 
estimates may be very heterogeneous; 2) Tx was not avail- 
able for the poorer clusters (N 3a z s < 20); and 3) the small 
sample size used in that work makes the results sensitive to 
the exclusion of a few points. 

We repeated the analysis presented in Figure 11 but 
only for clusters with N ff(l ; s ^ 20 and found that the scatter 
of the relations decreased only by 5%. However, the slope 



exhibited a significant variation (w 27%). So , the strong 
segreg ation present in the Tx-N 9 „| s relation of lLopes et al.l 
( 2006) may be explained by the small sa mple size and the ex - 
clusion of poor systems (N 3a ; s < 20). In lLopes et al.l {2009), 
we checked that the Tx values derived from BAX represent 
a consistent data set, which implies that the main source of 
slope variation detected is the absence of low mass systems. 



3.5 Comparison of scaling relations derived with 
optical and X-ray mass estimates 

As described in lLopes et al.1 (|2009l ) we have searched 
BAX (Base de Donnees Amas de Galaxies X, 
|http:/ /bax.ast. obs-mip.fr/ ) for counterparts to the 183 
clusters used in this work. The search was restricted to 
objects at z < 0.12 with X-ray temperature measures. We 
found 282 clusters in BAX, of which 21 are common to 
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Figure 6. Analogous to the previous figure, but showing only the 
relation M 2 oo-L5 200 for the 53 CIRS clusters. 



Figure 8. Analogous to the previous figure ("frame" back- 
ground), but showing only the CIRS clusters. 
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Figure 7. Analogous to Figure 5, but showing only the relation 
M200-Lx 200 f° r the full sampe and with the "frame" background. 



our sample. For these clusters, we employed the M 2 oo-Tjf 
relation given by equation 3 of POP05 to estimate masses. 
Our goal is to check if the mass calibration gives the same 
results when using either the optical or X-ray masses. Of 
the 183 NoSOCS plus CIRS systems, temperatures are 
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Figure 9. The relations ap-hf^ (top panels) and M 2 oo-Lo°? 
(bottom panels) are shown for all clusters (left) and systems with- 
out significant evidence for substructure (right). We apply the 
Dressler-Schectman (or A) substructure test to the spectroscopic 
sample of 170 clusters with at least 5 galaxies (with spectra) 
within R200- Of these, 39 (~ 23%) clusters shown strong signs 
of substructure. 
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Figure 10. Analogous to the previous figure, except for the data 
set and substructure test employed. Now we use the symmetry (or 
0) substructure test applied to the photometric data (m* — 1 
m r ^m* + 1 ) for the 179 clusters with at least 5 galaxies within 
R200- This is a two-dimensional (2D) test. Out of the 179 clusters, 
43 (~ 24%) systems are found to have substructure. 



generally available only for the more massive clusters, with 
g p > 4 00 km/s or M 2 oo > 10 14 M©. 

In I Lopes et all (|2009l ) we show that the temperature 
values selected from BAX represent a consistent data set. 
There, we compare the values listed in BAX to the ones 
available in RD06, finding that most systems agree within 
10%. We also check that the masses obtained from the dy- 
namical analysis of the optical data generally agree within 
40% of the estimates derived from the M200-T.Y relation. 
The exceptions, most of times, are clusters affected by sub- 
structure. When using the optical masses derived with the 
caustic technique (RD06) the agreement is not as good, with 
the optical masses being lower than the X-ray values in the 
low-mass regime. 

In this section, we check if the scaling relations obtained 
in this work are similar to the ones using mass derived from 
the tight connection to Tx ■ As the velocity dispersions and 
masses obtained with the caustic technique are provided by 
RD06 we also compare our results to those based on their 
estimates. Figures 12 and 13 show the M2oo-Lopt and M200- 
Lx relations, respectively. Lx is the value derived in this 
work with the "frame" background, using RASS data. Both 
figures consider M200 as determined from Tx (equation 3 
of POP05). These results are summarized in the first four 
lines of Table 7, where we also list the fit parameters for 
the M2oo-N sa ; s relation and for M200-LX (considering the 
"annulus" background) . 

From inspection of this table, and considering the range 
of clusters sampled, the results are in very good agreement to 
those based on the virial mass and restricted to the richer 
clusters. In Tables 1, 3, 4 and 5 the CIRS sample repre- 




Figure 11. Analogous to the previous figure, except for the sam- 
ple of galaxies considered. The /3 test is used again. However, we 
consider the values of L op t computed within 0.5 h" 1 Mpc and the 
substructure results obtained within 1.5 h — 1 Mpc. The choice of 
optimal sizes of the fixed aperture for counting gala xies and the 
way w e estimate substructure are both explained in lLopes et al.l 
jjroj ). The f3 test is applied to all 183 clusters of the NoSOCS 
and CIRS sample. We find that 64 (~ 35%) clusters show signs 
of substructure. 
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Figure 12. The relation between mass estimated from X-ray data 
and optical luminosity. 
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Figure 13. The relation between mass estimated from Tx and X- 
ray luminosity. Lx is the one determined in this work, considering 
the "frame" background. 



sents the richer clusters. Although the sample with masses 
derived from Tx (Table 7) is less than half of our rich sam- 
ple, we conclude that using the optical or X-ray mass leads 
essentially to the same scaling relations. In particular, the 
slopes are always concordant within 1-a. It is imperative to 
compare the results in Table 7 to those of the rich samples 
of the previous tables, which have approximately the same 
mass range. 

Note also that the scatters of the relations in Table 
7 are generally smaller than what we found with the virial 
masses. However, the results agree within 1-a. The improve- 
ment is more pronounced only for the relations involving 
Lx, reaching half the scatter for the "frame" background 
(but still within 2-er). These more accurate results are due 
to the smaller sample size of the relations based on the mass 
obtained from Tx (only 21 clusters). The most important 



conclusion from this exercise is that N 



gals 



-'opt 



and Lx can 



be used for mass calibration and the results are independent 
to the way mass is estimated, from optical (virial analysis) 
or X-rays (using the M-Tx relation). 

We decided to investigate if these conclusions hold for 
the parameters derived from the caustic technique (RD06). 



Figures 14, 15, 16, and 17 show the M200-N 



gals i 



M, 



-'opt , 



M200-LX and ap-Lopt relations, considering mass and ap 
obtained with the caustic technique. These parameters are 
provided for all CIRS clusters in RD06. Richness, optical and 
X-ray luminosities are estimated in the present work. These 
relations, as well as a\p-N sa ; s and ap-hx, are summarized 
in the last eight rows of Table 7. 

When comparing these relations to the results in Tables 
1, 3, 4 and 5 (only for the CIRS systems) the normalizations 
show better agreement than above. This is due to the fact 
that all the comparisons use only the CIRS sample. The 



slopes agree within 2-a for the relations involving ap and 
within 1-a for those regarding M2oo- However, this only hap- 
pens because the uncertainties in the fit parameters of the 
relations involving mass are much larger than the ones re- 
garding ap. The uncertainty in the slopes of the relations 
using the caustic mass is also larger than what we found 
with the virial mass, by as much as a factor of three. So, 
although the relations based on our masses (or the X-ray 
derived masses) agree within the errors to those obtained 
with the caustic method, the nominal values are very differ- 
ent. When using the caustic parameters, the relations about 
ap have smaller slopes compared to ours. The opposite is 
true for the relations regarding M2oo- 

We also find that the relations based on the caustic 
values are noisier (with larger scatter) than ours, especially 
for results involving mass. This is easily seen in Figures 14-17 
and Table 7. For instance, for the M2oo-L^, 2 t 00 relation the 
scatter in mass at fixed luminosity is (0.74 ± 0.18) when 
using the CIRS values, but only (0.36 ± 0.06) for the virial 
masses determined in the current work. 

Such findings are explained by the fact that our veloc- 
ity dispersions are in good agreement with those of RD06, 
as shown in paper I. However, the same is not true when 
comparing masses. Those masses estimated from Tx also 
disagree with the caustic values. RD06 claim concordance 
between the caustic results and their virial masses (as well 
as the X-ray estimates). However, comparing our masses and 
the caustic values (RD06) to the X-ray estimates, we find 
a better correlation with our values. Therefore, we conclude 
that our estimates better represent the cluster potentials, 
as discussed in Section 6 of paper I. There, we also argue 
that the interloper removal procedure should not be respon- 
sible for the different results (as the velocity dispersions are 
similar between our work and RD06). 

The culprit for the different mass estimates lies in the 
discrepant values of R200, which in RD06 are derived from 
the caustic mass profile. In paper I we estimate R200 as 
a by-product of the virial analysis an d the assumption of 
an NFW profile l|Navarro et alii 19971 ). This profile is well 
matched to the caustic mass profile for only half of the CIRS 
sample (RD06) , which may explain the different results and 
the larger scatter for the caustic based scaling relations. 
Note also that RD06 use R200 determined by the caustic to 
compute the virial masses. That choice helps improving the 
agreement they find between the caustic and virial masses . 
The conclusions we reached in paper I (|Lopes et al.l l2009h 
are corroborated by the scaling relations shown above. 



3.6 The mass-to-light ratio 

In this section, we present the mass-to-light ratio (M/L) 
of clusters of galaxies. That has been used before for es- 
timating the mass density of the universe. The typical 
value of M/L for rich clust ers is M/L ~ 300 h M Q /L 
l|Bahcall fc Comerfordl 120021 ). Previous works have shown 
that M/L increases with the size of the bound systems, from 
galaxies to clusters. The dependence of the cluster M/L ra- 
tio on cluster mass can be par ametrized by a relation like 
M/L oc M°, with a ~ 0.20 jBahcall fc Comerfordl [20021 ; 
iPopesso etld1l2005l ). 

In Figure 18, we show the mass-to-light ratio obtained 
in this work. We recomputed L op t considering the interval 
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Figure 14. The connection between the caustic mass (RD06, 
estimated from the optical data) and richness. The relation is 
exhibited for the 53 CIRS clusters. 
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Figure 15. Analogous to previous figure, but showing L op t. 

of m*-5 to m*+5. Using the values obtained in paper I and 
considering the same range used for the richness estimation 
(m*-l to m*+2), would lead to M/L values that are too 
high. We do not consider an incompleteness correction to 
the luminosity function (LF) for galaxies fainter than m*+5 
as the correction is around 5% and uncertainties in the slope 
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Figure 16. The relation between the caustic mass (RD06) and 
X-ray luminosity. Lx is obtained with the "annulus" background. 

at the faint end of the LF could bias the L op t values. So, 
integrating the LF to infinity does not significantly change 
L opt . The range we considered is already deep enough (being 
five magnitudes fainter than m*). In Figure 18, the solid line 
shows the relation listed i n Table 3 for ALL cl usters, while 
the other two are based in IPopesso et "all (|2005h . The dotted 
line considers their "optical" sample (69 clusters with optical 
masses), while the dashed line is for their "enlarged" sample, 
with 102 clusters comprising the optical sample plus clusters 
with masses estimated from the M-Tx relation. All the lines 
are normalized to an M/L ratio of 200. 

As we can see from Figure 18 our results show good 
consistency to previous findings, reinforcing the dependence 
of the M/L ratio to the cluster scale. Note that the power 
of the relatio n we find is a = 0.1 0, smaller than the value 
of 0.19 from IPopesso et al.l (|2005h ("optical" sample; dot- 
ted line). However, both results are consistent within 1-a 
and we find nearly the same slope as represented by their 
"enlarged" sample (a — 0.09; dashed line). The mean M/L 
ratio we find is ~ 210, in excellent agreement to previous 
studies l|Bahcall fc Comerfordl l2002f l. In a future work we 
plan to investigate the LF function of clusters in detail and 
in conjunction to the mass-to-light ratio and its dependence 
on cluster mass. 



4 COMPARISON TO THE LITERATURE AND 
DISCUSSION 

The results in Tables 1, 3, 4 and 5 indicate a good agree- 
ment with most findings in the literature. In Table 8 we 
summarize our main results, considering all clusters, or only 
the CIRS systems, as well as the findings of other authors. 
In particular, there is good consistency between the slopes 
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Figure 17. The relation between velocity dispersion estimated 
with the caustic technique and optical luminosity. The relation is 
shown for the 53 CIRS clusters. 




Figure 18. The mass-to-light ratio for the 127 NoSOCS plus 
53 CIRS clusters. The solid line shows the relation listed in Ta- 
ble 3 for ALL clust ers. The other two results shown are from 
iPopesso et al. | d2005h . The dotted line considers their "optical" 
sample, while the dashed line is for their "enlarged" sample. All 
the lines are normalized to a M/L ratio of 200. 



and scatter of POP05 (also based on SDSS and RASS) with 
our results, in the worst cases within 1.5-ct. RD06 claim to 
find good consistency with POP05. However, as seen in the 
previous section that is at the cost of larger errors for the 
fit parameters and much noisier relations. Note that our 
results based solely on CIRS (the most massive systems) 
always have smaller slopes (an issue we discuss furth er be- 
low) . Using an enlarged data set IPopesso et al.l l|2007l ) found 
relations consistent to POP05. 

The works of POP05 and RD06 are the easiest compari- 
son to ours, since they are also based in SDSS data. However, 
other results, b ased on different surveys and w avelengths can 
be informative. iLin. Mohr fc Stanford! (|2004 ) found M200 oc 
L x?H200 and M200 oc Ng a 1 ; 5 fl2O0 , using Two Micron All Sky 
Survey (2MASS) data. Our slopes agree within 1.5 and 1-er, 
respectively. Note that we consider their results after exclud- 
ing the BCG, as the y claim better agreemen t on the powers 
of the two relations. iKochanek et al.l ((2003) found M200 oc 
N»'|g 6 , also using the 2MASS K-band, a result that agrees 
with ours within 1.5-q. These fi ndings corroborate the con- 
clusion of IPopesso et al] 1120071 ). who reports no significant 
difference among results in different SDSS bands. So, even 
considering the 2MASS K-band it seems that the results are 
similar. 

On what regards other works using Ljy as a mas s 
tracer, we find th a t the slopes of I Popesso et al.l d2005l); 
IStanek et all l|2006l ); lvikhlinin et "all (|2008l 1 and lEttori etaD 
(HI) are consistent to ours. The scatter is consistent 
(within 1-a) to th e works of IPopesso et al. (20051) and 
IStanek et all (120061). whi le the scatter from IVikhlinin et alj 
i|200Sl ) and lEttori et all (|2004l ) are smaller than ours (con- 
sistent only within 3- c r). Note that the fir st two works above 
l|Popesso et al.1 120051 ; IStanek et all 120061 ) are based in the 
RASS data, while the last two use the Chandra satellite. 
We find a very good agreement (within 1-a ) between the 
norma lization we measure and the results of IPopesso et al.l 
(2005), for the M-Lx relation. Differences in the normaliza- 
tion to the other works above might be due to the differ- 
ences in the cluster samples used, which may span different 
mass/luminosity regimes. The comparison of the different 
selection functions is very hard and is not our goal in this 
work. 

The small differences with POP05 still deserve a few 
comments. The scatter in M500 at fixed L^ t 00 is similar to 
ours (consistency within 1.5-cr). However, a closer look at the 
upper left panel of our Figure 4 and their Figure 5 shows 
that we have many more clusters at M500 <1O 14 M0 than 
they do. So, their results are based on a data set showing 
a cloud of points at M500 >10 14 M© and a few scattered 
systems below this value. Nonetheless, that small data set 
in the low mass regime is sufficient to put their results and 
ours at the same level. One odd aspect of their plot is the 
diffuse distribution of the high mass systems, showing high 
scatter in this mass regime. Our relations always show a 
reduced scatter towards higher mass. Their poor sampling 
at low masses can also be seen in Figures 11 and 13 of their 
work, where we note a deficit of clusters below ap = 400 
km/s. 

Considering this sampling issue in POP05 it is natural 
to compare their results with the ones we obtained using 
only CIRS (which also has mostly clusters with ap > 400 
km/s). Figure 3 shows the relation between ap and L^J 00 . 
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It is remarkable how low the scatter is in this figure. How- 
ever, the most interesting point is that the relations are 
slightly flattened for the CIRS clusters, when compared only 
to NoSOCS or all systems. So, it seems that the scaling re- 
lations have different behaviors in different mass regimes. 
They look steeper for the low-mass systems, becoming flat- 
ter for the high-mass clusters. Using CIRS we note that the 
agreement to POP05 is as good as before, despite the fact 
the slope is now smaller. 

iPopesso et al.l 120078) created a larger sample by adding 
130 Abell clusters with confirmed 3D overdensities in SDSS. 
Although this sample is larger than ours (217 clusters com- 
pared to 180) it still seems that there are few systems below 
<10 14 Mq (see their Figure 9). Most confirmed Abell over- 
densities are rich systems, which explain the figure (note 
that the open symbols are the X-ray systems, which con- 
tain the lower mass objects). The main conclusions drawn 
from these comparisons are: (i) the slopes of the scaling re- 
lations become flatter for high-mass samples, but the dif- 
fere nces are within 1-cr; (ii) although the samples of POP05 
and lPopesso et all (|20Q7l ) show fewer systems with low mass, 
their results are consistent with ours. That is due to the in- 
clusion of a few groups in their sample and to the small 
variation in the slope with mass. These results indicate that 
sample incompleteness has little effect on the scaling rela- 
tions. 

For the relations involving Lx we find the same trend 
as above, with the high mass systems showing flatter rela- 
tions. However, some clusters were not considered in the fits 
due to having only upper limits in the X-ray luminosity esti- 
mates. This eliminates points that generally have Lx < 10 43 
ergs/s in Figure 5 (lower right panel). These points would 
be located above the derived fit and would thus make the 
relations a little flatter. Our results with all clusters could 
be biased due to missing some low mass systems with Lx ~ 
10 43 ergs/s, making the relations steeper than it would re- 
ally be. However, considering the conclusions reached above, 
we might expect that incompleteness (especially in the low 
mass regime) has a minor impact in the relations. We also 
note that the results within R500 generally are the ones with 
smallest scatter when considering Lx- 

A result that deserves some comments is the third test 
shown in Table 2 (rows 5 and 6) where we consider all clus- 
ters at z ^ 0.25, considering that at z > 0.10 the spectro- 
scopic survey of SDSS is no longer complete. The main ef- 
fect exhibited is the flattening of the scaling relations. That 
happens because N ga is and L opt are still well determined 
(as they rely on the photometric data), but ap and mass 
are not (see discussion in paper I). If the percentage of clus- 
ters without enough spectra is large, the impact in the scal- 
ing relations may be important. So, the safest approach is 
to consider only objects that have complete spectroscopic 
sampling (as done here and by RD06). 

Another importa nt issue concerns t he interloper re- 
moval procedure ( se e | Woitak et al.l l2007h . As discussed in 
§4 of lLopes et al.l (|2009T 1 we performed several tests to 
optimize the "shifting-gapper" technique. The final code 
seems to work well for systems of all masses, as can b e 
seen from the phase-space diag rams llLopes et alj I2009T) . 
PO P05 considers the method of iKatgert et all l|2004h Tsee 
also lden Hartog fe Katgertlll996l ). which also combines the 
position and velocity information, but makes assumptions 



regarding the dynamical state of the cluster. RD06 use the 
caustic technique to determine the mass profile of the clus- 
ters, which involves removal of int erlopers in underdens e re- 
gions of the phase-space diagram. IWoitak et all l|2007h ap- 
plied several different techniques for interloper removal to 
real and simulated data. They conclude that some meth- 
ods have poor performance compared to others, but none 
of the methods cited above are among the poor performers. 
Whatever the approach, it must be sufficiently restrictive 
but not overly lax, otherwise many cluster members will be 
discarded or few interlopers will be rejected, leading to bi- 
ased values of a p and mass. While there is no particular rea- 
son to chose one method over another, the procedure should 
work for a wide variety of systems (as shown for the method 
adopted here). We based our choice on the simplicity of the 
method and on the avoi dance of assumptio ns regarding the 
dynamics of the cluster (IFadda et al.lll996f ). 

Results from a test with a less rigorous criterion (elimi- 
nating fewer interlopers) are shown in rows 3 and 4 of Table 
2. The comparison of Tables 1 and 2 reveals no systematic ef- 
fect on the intercept and slope of the scaling relations when 
relaxing the interloper removal procedure. The scatter of 
the relations increases slightly, but the results are still con- 
sistent with our original relations and to those in the litera- 
ture. So, we conclude that our method is robust, rendering 
similar lists of cluster members than other methods. This is 
reflected in the good agreement with the velocity dispersions 
of RD06 and the comparison to the scaling relations of other 
authors jKochanek et aill2003l ; iLin. Mohr fc Stanfordll2004l ; 
IPopesso et al.ll2005h . 

The Halo Occupation Distribution (HOD) relates the 
number of galaxies within the virial radius and the associ- 
ated mass, expressed as N ga i,R200 oc M^oo- Hierarchical mod- 
els of structure formation predict that the number of subha- 
los within a system is directly proportional to the mass of the 
parent halo. In other words, a = 1. However, different mech- 
anisms (such as a decreasing efficiency of star formation, or 
an increased merger rate) would imply a decreasing number 
of galaxies per given mass i n higher mass halos , agreeing 
with a < 1 (see discussion in IPopesso et al.ll2007T ). Here, we 
do not intend to investigate the HOD. Our main goals are to 
provide calibrators to the cluster mass and compare those to 
determine which one traces mass more accurately. However, 
it is worth mentioning that the connection between mass 
and richness we found is indeed consistent with a < 1, as 
most recent results in the literature. In a future work we 
plan to study in detail the HOD and the factors that could 
make a < 1. 

Finally, we would like to comment that when compar- 
ing the performance of N s(l ; s , L opt and Lx as mass tracers, 
we find that the optical parameters show relations slightly 
tighter than L X - The scatter of M 500 at fixed N^f, 00 , L^f t 00 
and Lf 500 is (0.43 ± 0.03), (0.44 ± 0.04) and (0.47 ± 0.05), 



respectively. The scatter of M200 at fixed N^ 2 , 00 , L^ 200 



and 

L x 2m ' is (0.46 ± 0.04), (0.49 ± 0.04) and (0.56 ± 0.06), 
respectively. In all cases, the scatter and the slope of the 
relations are reduced for high mass systems. POP05 argued 
that the poorer performance of Lx is likely due to the vari- 
ation in the compactness of the galaxy clusters. They also 
investigated if the scatter in the M-Lx relation could be due 
to cooling core effects, but they find that this accounts for at 
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Figure 19. Correlation between X-ray luminosity and richness, 
both measured within R200. Lx is estimated with the "frame" 
background. 

most 3% of the scatter. The low signal to noise in RASS (es- 
pecially for the low-mass systems) can also contribute to the 
higher scatter found in the relations involving Lx ■ However, 
it is important to note that here we quote the errors in the 
scatter, finding that they are all consistent (in the worst case 
within 1.5-cr) for the three mass proxies considered (N ga i s , 
L opt and Lx). 



5 CORRELATION BETWEEN OPTICAL AND 
X-RAY PROPERTIES 

The interplay between global properties obtained in the op- 
tical and X-ray regimes are crucial for understanding the 
complex physics present in galaxy clusters. In recent years 
a number of works have been dedicated to the comparison 
of X-ray and o ptical catalogs or the construction of com- 
bined samples jDonahue et alJ l200ll. |2002|; Gilb ank et al.l 
l2004l : |Popesso et al.ll2004l . l2005l : lLopes et alJ2006f ). The abil- 
ity to predict the X-ray luminosity or temperature from op- 
tical parameters, and vice-versa, is also important for fu- 
ture surveys, conducted only in one regime. As a by-product 
these comparisons may also reveal unusual clusters (such as 
X-ray underluminous ones), which are very interesting for 
follow-up studies. 

In this section, we investigate the correlations between 
optical (N ga i s and L op t) and X-ray (Lx and Tx) quanti- 
ties. Note that the X-ray luminosity is estimated through 
an iterative proce dure that assumes the Lx-Tx relation of 
iMarkevitchl (|l998h (see paper I for more details). That re- 
lation is close to the self-similar expectations (without any 
correction for cooling flows in Lx or Tx)- This is the same 
relation adopted to estimate Lx for NORAS, in which Tx oc 
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Figure 20. Correlation between X-ray and optical luminosities, 
both measured within R200. Lx is estimated with the "annulus" 
background. 
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L^ ■ Figures 19 and 20 show the connection of Lx to N ga ( 3 
and L op i, respectively. In the first relation Lx is the value 
obtained with the "frame" background, while the "annulus" 
background is used in the latter. In both plots, all mea- 
surements are performed within R2oo- Figure 21 shows the 
correlation between Tx (given by BAX) and L^ p 2 t 00 for 21 
clusters. 

Table 9 summarizes the fitting parameters for the rela- 
tions between X-ray and optical measurements within 0.50 
h _1 Mpc, R200, and Rsoo- Note that if one does not have a 
measure of a radius that scales with mass (R500 or R200), 
0.50 h _1 Mpc represents the optimal fi xed aperture for com - 
pari ng optical and X-r ay properties. iPopesso et al.l (|2004 ) 
and lLopes et alJ (|2006h tested several different fixed aper- 
tures and found the scatter to be the smallest when using 
0.50 h" 1 Mpc. The results in Table 9 are listed for the two X- 
ray background types considered ("annulus" and "frame"). 
In Table 10, we list the results involving Tx (21 clusters 
found in BAX) for all three apertures used for computing 
N ga!s and L opt . 

Our results are in good agreement with others in the lit- 
erature. For instance, POP05, also using data in the SDSS 
r-band, found similar results to ours: Lx,H200 oc Ll' p 7 ^ R200 
and Tx oc Lo^^2oo- The results within R500 and with a 
fixed metric (in this case for the i-band) are also close to 
ours. In the worst ca ses consistency is found within 1.5-cr. 
iDonahue et~aH l|200lD found similar relations but using a 
different richness definition, which hampers a direct com- 
parison to our findings. 

We also find a good agreement t o the Tx-N 9a ; s and 
Tx-L pi relations of lLopes et al.l |2006;). The differences are 
within 1-ct. The relation s invo lving Lx are not as similar. 
Note that in lLopes et all (|2006l ) we also force a richness cut, 
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Figure 21. Correlation between temperature and optical lumi- 
nosity (measured within R200). 



Figure 22. Correlation between radius (R200) and N^ 2 ^ . 



with N aa i a ^ 10 (richness within 0.5 h _1 Mpc) . lLopes et ail 
(2006) also considered optical estimates from DPOSS and X- 
ray values from BAX. The heterogeneous nature of the lumi- 
nosities derived from BAX could contribute to the discrep- 
ancies. However, the relations involving Lx are much steeper 
and thus more sensitive to the sample used. If we only con- 
sider the most rich clusters the agreement is much better 
than with the full sample. We also notice that the correla- 
tions fcamcHn thepresent work are in excellent agreement 
with lGal et ail (|2009| ). who used DPOSS data for measuring 
and Lopt within 0.50 h _1 Mpc and directly estimated 



N 



gals 



Lx from RASS, as done here. 

If the density profiles for dark matter and intra-cluster 
gas are self-similar, the following relations hold: M oc T 3//2 , 
Assuming mass traces light (con- 



T oc \}Jr 2 and M oc L x 



3/4 



stant M/L j,i) we would expect that L opt oc L 3 ^ 4 and L op t oc 
T 3//2 . If there is also a strict proportionality between L opt 
and richness, the same relations above are valid for N ga z 3 . 
Our results are consistent within 2-cr with these relations. 
However, if we assume that T oc L 2 / 5 (the typically observed 
result, iDavid et al.l|l993t ) we would have L opt oc L 3 ^ 5 , con- 
sistent with our findings at 1-a. 



parent halo, namely a = 1. However, different mechanisms 
such as a decreasing efficiency of star formation, or an in- 
creased merger, and destruction rate of galaxies, would im- 
ply a decreasing number of galax ies for higher mass halos, 
suggesting a < 1 (see discussion in lPopesso et al-lfeoOTl) . Fig- 
ure 22 shows the connection between R200 and N^ 2 ^ , while 
Table 11 has the solutions for the fits regarding the R200- 



N 



0.5Mpc/h 
gals 



R200-N 



gals , R 500-N gais 



0.5Mpc/h 



and Rsoo-N^ 



We find R200 oc N l ' a J 



0.39 

gals, fl200i 



which indicates that a 



0.86. The connection between N 



H200 
gals 



and M200 shown in 



Table 1 indicates a = 0.93, agreeing within 1-a to the 
finding above. The connection between R200 and N ga i s 
we found is al so very close to the result obtained by 
iGal et ail (120091) who est imated R200 photometrically (simi- 
lar to Hansen et al .120051 ), but counted galaxies as in the cur- 
rent work. This agreement corroborates the findin gs and re- 
liabili ty of the photometric R200 estimates given in Gal et ail 
(|2009h . 

In Table 11 we also list the results relative to optical lu- 
minosity. We see that the relations for L opt are consistent to 
the ones based on richness. These results are also useful for 
compa rison to the slopes of the relations from lPopesso et"ail 
(|2005h . In the SDSS r-band they found that R500 oc L° p 4 t ° fl500 
and R200 oc L° p 4 * H2 00! m good agreement with our findings. 



6 CONNECTION BETWEEN RICHNESS AND 
RADIUS 

We know that the radius and mass of a cluster scale as 

1/3 

R200 oc M 2 q and the number of galaxies within R is linked 
to mass as ~N ga is,R200 oc Mfoo- If N 9a ; s is measured within a 
fixed aperture the n N oai3 oc , p ^ n „ , with /3 ~ 0.50 — 0.65 
(see discussion in lHansen et al.lT2005h . Hierarchical models 
of structure formation predict that the number of sub-halos 
within a system is directly proportional to the mass of the 



7 CONCLUSIONS 

We have used a sample of 127 NoSOCS plus 56 CIRS 
galaxy clusters to investigate scaling relations at low red- 
shift (z ^ 0.10). For every cluster we previously determined 
(in paper I) the velocity dispersion {up), physical radii (R500 
and R200), masses (M500 and M200), richness (Ng a!s ), opti- 
cal and X-ray luminosities (L opt and Lx). The last three 
parameters are estimated within 0.5 h _1 Mpc, R500 and 
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R.200- Substructure estimates are also available for nearly 
all clusters. We estimated the presence of substructure from 
the galaxy distribution in two and three dimensions, using 
the f3 and A tests, respectively. For the CIRS systems we 
also consider the values of ap and M200, independently de- 
termined with the caustic method (RD06). For a subset of 
21 clusters we have Tx values from the literature and es- 
timated their masses using a M200-TX relation. The main 
conclusions we reach are: 

(i) Richness, optical and X-ray luminosities correlate well 
with cr p and mass. However, the results indicate that the 
slope and scatter of the relations are higher when using a 
fixed aperture to compute N ga i s and L op t- So, the most ac- 
curate relations are achieved when considering a physical 
radius, such as R500 and R2<x>- In general, the most accurate 
results are obtained within Rsoo- 

(ii) The scaling relations derived only with the CIRS sam- 
ple are flatter than the global relations. This indicates that 
the most massive systems scale differently than the poorer 
ones. However, the results for the poor and rich systems are 
still consistent within 1-a. 

(iii) The scaling relations show no significant modifica- 
tion when considering only clusters without substructure. 
The new relations, as well as the scatter, are consistent 
with the original ones. This result holds if we employ a two- 
dimensional (/3) or a three-dimensional (A) test. In partic- 
ular, the 2D test leads to the same scaling relations when 
we use the spectroscopically or photometric selected galaxy 
samples. 

(iv) The comparison of the scaling relations obtained with 
optical (virial) and X-ray masses (from the M200-TX rela- 
tion) point to very similar results. In other words, the mass 
calibration with mass estimates from different wavelengths 
are equivalent. This is one of the main results of this work. 

(v) As the velocity dispersions and masses computed with 
the caustic technique are available in RD06 we have also de- 
rived the scaling relations with these parameters (consider- 
ing the values of N g(l ; s , L op t and Lx computed here). Our 
findings are in line with the conclusions of paper I. Although 
the slope of the relations agree within 1-a that is only true 
because the fits obtained with the caustic parameters have 
a very large scatter. For the IVboo-L^, 200 relation the scatter 
in mass at fixed luminosity is (0.74 ± 0.18) when using the 
CIRS results, and only (0.36 ± 0.06) for the virial masses 
obtained in the current work. The situation is not as critical 
for the relations based on ap, indicating that the interloper 
removal procedure is not the answer for the observed dis- 
crepancies. Instead, these are due to the procedure used for 
the mass estimation, derived from the caustic mass profile. 
Uncertainties in this profile lead to different values of R200 
and as consequence, biased results for M200 (see discussion 
in paper I). 

(vi) We find a good agreement with most of the results 
in the literature, even those derived from other bands (such 
as 2MASS K). So, the comparison between optical and X- 
ray properties, the connection between R200 and N go ; s , and 
the mass calibration performed with different cluster proper- 
ties, all agree well with previous findings. This als o indicates 
the in terloper removal procedure we employed (|Lopes et al.l 
120091 ) is robust, as other authors we compared to employ 



different techniques for selecting cluster members, or even 
determine mass from other wavelengths. 

(vii) The scaling relations based on clusters at 2 < 0.25 
arc flatter than our original results (considering only objects 
at z < 0.10). However, the results of POP05 — who use 
clusters in SDSS at 2 < 0.25 — are consistent with our find- 
ings. We argue that the fraction of higher-2 systems in their 
sample is probably small, contributing little to their results. 
However, we acknowledge the relevance of this issue and con- 
clude that clusters' velocity dispersion and masses must be 
estimated from complete spectroscopic samples (reaching at 
least M* + 1; see paper I). Using a large sample of clusters 
having spectra for only the bright members can severely bias 
ap and mass, thus affecting the scaling relations. 

(viii) The main result of this paper regards the mass cal- 
ibration of galaxy clusters. We show that richness, L op t and 
Lx can reliably be used for mass estimation in the nearby 
universe. Particularly, the optical properties provide slightly 
more accurate relations in the present work, which is prob- 
ably due to the use of RASS in the X-ray regime (shal- 
lower than the optical data from SDSS). When considering 
all clusters, the scatter of M 50 o at fixed Nf a 5 , 00 , L^ p °° and 
Lf 500 is (0.43 ± 0.03), (0.44 ± 0.04) and (0.47 ± 0.05), re- 
spectively. The scatter of M200 at fixed Nf a 2 , 00 , L^ p 2 t 00 and 
Lf 200 is (0.46 ± 0.04), (0.49 ± 0.04) and (0.56 ± 0.06), re- 
spectively. For the richer (CIRS) systems we found the scat- 
ter of M500 at fixed Nfi) 00 , L^°° and Lf 500 is (0.33 ± 0.05), 
(0.38 ± 0.05) and (0.48 ± 0.06), respectively. These findings 
are in accord with those of POP05, who found that L op t 
is a slightly better mass tracer than Lx (also considering 
SDSS and RASS). This is a very important conclusion, be- 
cause it tells us that with accurate single band photometry 
we can reliably estimate the mass of galaxy clusters, a key 
result for studying the cluster mass function at low redshift. 
Our work also indicates that the spectroscopic follow-up of 
a few dozens of clusters to high redshifts (2 ~ 1) can be 
used for the understanding of the evolution of the scaling 
relations and to trace mass at high— 2. This is a crucial step 
for self-calibration methods aiming to constrain the dark en- 
ergy from the evolution of the cluster population with cosmic 
time. 

Note that we do not state that N ga i s or L opt are better 
mass tracers than Lx- We only say that based on RASS 
and SDSS that is the case. That is due to the fact that 
RASS is shallow compared to SDSS. A new X-ray survey, 
for instance based in eROSITA, would provide accurate X- 
ray luminosities for large samples of clusters. However, as 
we showed here, richness and optical luminosity are obser- 
vationally cheap parameters that can also work as a mass 
proxy. This is an important result for future large sky sur- 
veys such as DES, Pan-STARRS, LSST and UKIDSS. If the 
mass-calibration relation and its evolution are known, these 
surveys can provide accurate estimates of richness or L opt 
and thus the mass of clusters at high redshifts. 
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APPENDIX A: UPDATE VALUES OF X-RAY 
LUMINOSITIES 

We fixed a bug in the code used to estimate X-ray luminos- 
ity and its associated error. Hence we update these values 
in the table Al below, which replaces Table 5 from paper 
I. The meaning of all columns is as in paper I. In the first 
column we give the cluster name; then in the next three 
columns we list the values of Lx (and its associated error) 
for the three backgrounds ( "annulus" , "frame" and "boxes" , 
respectively). In the fifth column we show the X-ray tem- 
perature measure from BAX (when available), while the in- 
terpolated temperature obtained from the Lx-Tx relation 
is shown in the sixth column. The last column indicates 
whether we had a significant detection (SD), an upper limit 
(UL), or if the measurements should be taken with concern 
(XX). The cluster is marked as XX if it is too close to a 
border, or if the background error or net count rates could 
not be determined. The last two columns of the table are 
obtained with the "annulus" background. 



NoSOCS in SDSS. II 21 



Table 1. Relations between op, M500 and M200 to N ga ; s , using three different apertures, 0.50 
h _1 Mpc, R500 and R200- In all tables the linear fit is of the form ln(Y) = A + B X ln(X/C), 
where X and Y are listed in the first column. Masses are in units of 10 14 Mq and velocity 
dispersion in km s — 1 . The results are also shown for three different samples, using the 127 
NoSOCS clusters, fifty-three CIRS clusters, and the combination of the two samples with 180 
clusters (second column). The intercept (A) and slope (B) are shown in the third and fourth 
columns, respectively. The scatter in the Y parameter at fixed X is in the fifth column. The 
total number of clusters for the first fit is shown as Ntot and the number of clusters after 
a 3-ct clip is given by N use . The pivot point depends on the sample, being C = 25 for the 
NoSOCS and full samples and C = 60 for the CIRS sample. 



Relation 
X 


Y 


Sample 






A 








B 










iY\X 


N tot 


N^ se 


N 0.50 

gals 




NoSOCS 


5 


.88 


± 


0.03 





.68 


± 


0. 


04 





.29 


± 


0.02 


127 


122 


N-R500 
gals 


p 


NoSOCS 


5 


.80 


± 


0.03 





,54 


± 


0. 


03 





.23 


± 


0.02 


127 


124 


N-R200 
gals 


Op 


NoSOCS 


5 


.70 


± 


0.03 





.56 


± 


0. 


03 





26 


± 


0.03 


127 


125 


N 0.50 
gals 


M200 


NoSOCS 





75 


± 


0.09 


1 


.45 


± 


0. 


16 





69 


± 


0.06 


127 


124 


N-R500 
gals 


M500 


NoSOCS 





53 


± 


0.06 


1 


01 


± 


0. 


07 





.53 


± 


0.05 


127 


126 


N-R200 
gals 


M200 


NoSOCS 





.36 


± 


0.07 


1 


09 


± 


0. 


09 





.60 


± 


0.07 


127 


126 


N 0.50 
gals 


Op 


CIRS 


6 


.49 


± 


0.03 





.49 


± 


0. 


05 





.17 


± 


0.02 


53 


51 


N-R500 
gals 


Op 


CIRS 


6 


.29 


± 


0.02 


0. 


40 


± 


0. 


03 


0. 


.14 


± 


0.02 


53 


51 


N-R200 
gals 


Op 


CIRS 


6 


.20 


± 


0.03 





41 


± 


0. 


03 


0. 


.15 


± 


0.01 


53 


51 


N 0.50 
gals 


M200 


CIRS 


2 


.15 


± 


0.08 


1 


36 


± 


0. 


19 





46 


± 


0.08 


53 


52 


N-R500 
gals 


M500 


CIRS 


1 


,54 


± 


0.05 


0. 


.95 


± 


0. 


09 





.33 


± 


0.05 


53 


52 


N-R200 
gals 


M200 


CIRS 


1 


.39 


± 


0.07 


1 


00 


± 


0. 


11 





.37 


± 


0.05 


53 


52 


N 0.50 
gals 


op 


ALL 


5 


.91 


± 


0.02 





.69 


± 


0. 


03 





.25 


± 


0.02 


180 


171 


N-R500 
gals 


Op 


ALL 


5 


.81 


± 


0.02 





52 


± 


0. 


02 





.19 


± 


0.01 


180 


175 


N-R200 
gals 


op 


ALL 


5. 


.69 


± 


0.02 





53 


± 


0. 


02 





.22 


± 


0.02 


180 


178 


N 0.50 
gals 


M200 


ALL 





79 


± 


0.05 


1 


.50 


± 


0. 


08 





.57 


± 


0.05 


180 


174 


N-R500 
gals 


M500 


ALL 





.57 


± 


0.04 


1 


.04 


± 


0. 


04 





43 


± 


0.03 


180 


178 


N-R200 
gals 


M200 


ALL 





37 


± 


0.04 


1 


.07 


± 


0. 


04 





.46 


± 


0.04 


180 


177 



Table 2. Fit parameters for the op-N ga i 3 , M200-N ga i 3 relations within R200 for three more 
cases, all based in the NoSOCS sample (same pivot point as in Table 1). Masses are in units 
of 10 14 Mq and velocity dispersion in km s _1 . Rows 1 and 2 considers the original (instead 
of luminosity- weighted) coordinates. Rows 3 and 4 are for the case when we consider a more 
"relaxed" criteria for rejecting interlopers (described in § 4.1 of paper I), while rows 5 and 
6 have the results for all the 219 NoSOCS clusters at z ^ 0.25, instead of z ^ 0.10. All the 
columns have the same meaning as in Table 1. 



Relation 
X 


Y 


Sample 


A 


B 


a lnY\X 


Nt ot 




N-R200 
gals 


op 


NoSOCS 


5.72 ± 0.03 


0.54 ± 0.04 


0.29 ± 0.03 


128 


127 


N-R200 
gals 


M200 


NoSOCS 


0.38 ± 0.06 


1.05 ± 0.09 


0.60 ± 0.07 


128 


127 


N-R200 
gals 


Op 


NoSOCS 


5.77 ± 0.03 


0.58 ± 0.04 


0.29 ± 0.04 


130 


128 


N-R200 
gals 


M200 


NoSOCS 


0.51 ± 0.07 


1.16 ± 0.09 


0.66 ± 0.09 


130 


129 


N-R200 
gals 


Op 


NoSOCS 


5.63 ± 0.02 


0.43 ± 0.02 


0.26 ± 0.02 


219 


216 


N-R200 
oals 


M200 


NoSOCS 


0.21 ± 0.04 


0.83 ± 0.05 


0.61 ± 0.05 


219 


216 
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Table 3. Analogous to Table 1, but listing the fit parameters for the relations of ap, M500 
and M200 to L op t. Masses are in units of 10 14 Mq and velocity dispersion in km s _1 . As 
before, the pivot point depends on the sample, being C = 0.40 10 12 Lq for the NoSOCS and 
full samples and C = 1.10 10 12 L Q for the CIRS sample. 



Relation 
X 


Y 


Sample 






A 








B 








N tot 




T 0.50 
U opt 


Up 


NoSOCS 


5. 


.88 


± 


0.05 


0. 


.93 


± 





.10 


0.41 


± 


0.04 


127 


124 


t i?500 
'-'opt 


ap 


NoSOCS 


5 


.79 


± 


0.03 


0. 


.63 


± 





.04 


0.26 


± 


0.02 


127 


123 


t il200 
^opt 


a P 


NoSOCS 


5. 


.68 


± 


0.04 


0. 


,64 


± 


0. 


05 


0.38 


± 


0.05 


127 


125 


t 0.50 
^opt 


M200 


NoSOCS 





.76 


± 


0.13 


2 


.00 


± 


0. 


.32 


0.92 


-j- 


0.20 


127 


124 


t i?500 
^opt 


M500 


NoSOCS 


0. 


.51 


± 


0.06 


1 


.09 


± 





.07 


0.57 




0.06 


127 


125 


t i?200 
^opt 


M200 


NoSOCS 





.27 


± 


0.07 


1 


.18 


± 





.09 


0.63 


± 


0.08 


127 


123 


T 0.50 
^opt 


a P 


CIRS 


6 


.55 


± 


0.04 


0. 


.59 


± 


0. 


07 


0.20 


± 


0.03 


53 


52 


t i?500 
^opt 


a P 


CIRS 


6 


.31 


± 


0.02 


0. 


.41 


± 





.03 


0.14 


± 


0.01 


53 


50 


t i?200 
U opt 


a P 


CIRS 


6 


.22 


± 


0.03 


0. 


39 


± 


0. 


03 


0.15 


± 


0.02 


53 


51 


T 0.50 
^opt 


M200 


CIRS 


2 


.26 


± 


0.11 


1 


,48 


± 


0. 


22 


0.57 


± 


0.10 


53 


53 


t B500 
^opt 


M500 


CIRS 


1 


.57 


± 


0.05 





.99 


± 





.10 


0.38 


± 


0.05 


53 


53 


t -R200 
U opt 


M200 


CIRS 


1 


41 


± 


0.06 





98 


± 


0. 


.11 


0.36 


± 


0.06 


53 


51 


T 0.50 
^opt 


a P 


ALL 


5 


.85 


± 


0.03 


0. 


,78 


± 


0. 


05 


0.31 


± 


0.02 


180 


178 


t i?500 
^opt 


a P 


ALL 


5 


.76 


± 


0.02 





.56 


± 





02 


0.21 


± 


0.02 


180 


174 


t i?200 
^opt 


a P 


ALL 


5 


.65 


± 


0.03 


0. 


,55 


± 


0. 


02 


0.24 


± 


0.02 


180 


176 


T 0.50 
^opt 


M200 


ALL 





.68 


± 


0.07 


1 


70 


± 


0. 


,12 


0.72 


± 


0.07 


180 


178 


t i?500 
^opt 


M500 


ALL 





.49 


± 


0.04 


1 


.06 


± 


0. 


,04 


0.44 


± 


0.04 


180 


173 


T H200 
^opt 


M200 


ALL 


0. 


.28 




0.05 


1 


11 




0. 


,04 


0.49 


± 


0.04 


180 


174 



Table 4. Fit parameters for the ap-Lx, Msoo-Ljf, M200-LX relations using two different apertures, 
R500 and R200- The results consider the X-ray luminosity obtained with the "annulus" background 
(second column). Masses are in units of 10 14 Mq and velocity dispersion in km s _1 . First we show the 
results derived with all clusters, then only the CIRS systems are used. The remaining columns have the 
same meaning as in the previous tables. The pivot points are C = 0.20 10 44 erg/s for the full sampe and 
C = 0.55 10 44 crg/s for the CIRS sampe. 



Relation 
X 


Y 


Bkg Type 


Sample 






A 








B 








VlnYlX 


N iot 




t H500 


a P 


Annulus 


ALL 


5 


.90 


± 


0.03 


0, 


26 


± 





03 





,29 


± 0.04 


104 


104 


L fl200 


a P 


Annulus 


ALL 


5 


.85 


± 


0.04 


0, 


,27 


± 





.03 





,29 


± 0.03 


97 


97 


T -R500 


M500 


Annulus 


ALL 





68 


± 


0.06 





62 


± 





,04 





.47 


± 0.05 


104 


103 


t R.200 
^X 


M200 


Annulus 


ALL 


0, 


59 


± 


0.08 





67 


± 





06 





.56 


± 0.06 


97 


97 


Lfl500 


a P 


Annulus 


CIRS 


6 


21 


± 


0.04 





18 


± 





,04 





.26 


± 0.03 


52 


52 


L fl200 


a P 


Annulus 


CIRS 


6 


19 


± 


0.04 





19 


± 





,04 





.25 


± 0.03 


52 


52 


L H500 


M500 


Annulus 


CIRS 


1 


.39 


± 


0.08 


0, 


49 


± 





,08 





.48 


± 0.06 


52 


52 


L fl200 


M200 


Annulus 


CIRS 


1 


40 


± 


0.08 





52 


± 





08 





.49 


± 0.06 


52 


52 
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Table 5. Analogous to the previous table, but for the results based in the "frame" background. Masses 
arc in units of 10 14 Mq and velocity dispersion in km s _1 . 



Relation 
X 


Y 


Bkg Type 


Sample 






A 






B 




&lnY\X 


N iot 


N use. 


t H500 
h X 


a P 


Frame 


ALL 


5 


.92 


± 


0.03 


0.24 


± 


0.02 


0.28 ± 0.04 


98 


98 


t i?200 

L x 


a P 


Frame 


ALL 


5. 


90 


± 


0.04 


0.22 


± 


0.03 


0.28 ± 0.03 


90 


90 


T i?500 
L X 


M500 


Frame 


ALL 





.73 


± 


0.06 


0.59 


± 


0.04 


0.46 ± 0.04 


98 


96 


L fl200 


M200 


Frame 


ALL 





68 


± 


0.08 


0.61 


± 


0.06 


0.55 ± 0.05 


90 


89 


L fl500 


Up 


Frame 


cms 


6 


20 


± 


0.04 


0.18 


± 


0.04 


0.26 ± 0.03 


50 


50 


L fl200 


Up 


Frame 


cms 


6 


19 


± 


0.04 


0.18 


± 


0.03 


0.25 ± 0.03 


49 


49 


T i?500 
L X 


M500 


Frame 


cms 


1 


38 


± 


0.08 


0.49 


± 


0.07 


0.49 ± 0.06 


50 


50 


L fl200 


M200 


Frame 


cms 


1 


.40 


± 


0.08 


0.50 


± 


0.07 


0.50 ± 0.06 


49 


49 



Table 6. Fit parameters for the M2oo-N oa ; s and M200 _ L O pt for clusters with or without substructure. Masses 
are in units of 10 14 Mq and velocity dispersion in km s — 1 . First, we show (lines 1-6) the results considering 
an aperture of R200 f° r aJl clusters and for the substructure free systems. These are selected with the A or 
(3 tests, applied only to galaxies considered as cluster members (§4). Then we show the results considering 
the same aperture and the 2-D test but using the photometric sample (no restriction to galaxies with 
spectra). That is seen in rows 7-10. Finally, we list the results with a fixed metric of 1.5 h~ 1 Mpc for 
estimating substructure and using the f} test. Again, the photometric data are considered (rows 11-14). The 
parameters of each relation are listed in column 1, while the substructure test (when it is the case) is shown 
in column 2. The remaining columns have the same meaning as in the previous tables. The pivot points for 



all samples 


are C = 


25 when considering richness, and C = 


0.40 10 12 L Q 


for L op t. 






Relation 
X 


Y 


Sub-test 


Sample 


A 


B 


u lnY\X 


N tot 




N-R200 
gals 


M200 




ALL 


0.40 ± 0.05 


1.10 ± 0.05 


0.48 ± 0.04 


170 


167 


N-R200 
gals 


M200 


A 


Sub-free 


0.40 ± 0.05 


1.10 ± 0.06 


0.49 ± 0.05 


131 


129 


N-R200 

gals 


M200 


P 


Sub-free 


0.40 ± 0.05 


1.10 ± 0.05 


0.49 ± 0.04 


144 


141 


T -R200 
L opt 


M200 




ALL 


0.31 ± 0.05 


1.12 ± 0.05 


0.51 ± 0.04 


170 


163 


T -R200 


M200 


A 


Sub-free 


0.35 ± 0.06 


1.08 ± 0.07 


0.52 ± 0.06 


131 


127 


T K200 
^opt 


M200 


P 


Sub-free 


0.33 ± 0.05 


1.12 ± 0.05 


0.53 ± 0.04 


144 


138 


N-R200 
gals 


M200 




ALL 


0.38 ± 0.04 


1.09 ± 0.04 


0.48 ± 0.03 


179 


176 


jvj.R20() 
gals 


M200 


P 


Sub-free 


0.38 ± 0.05 


1.06 ± 0.05 


0.48 ± 0.04 


136 


133 


T H200 
^opt 


M200 




ALL 


0.28 ± 0.05 


1.12 ± 0.05 


0.52 ± 0.04 


179 


173 


T -R200 


M200 


P 


Sub- free 


0.33 ± 0.06 


1.07 ± 0.06 


0.51 ± 0.05 


136 


131 


N 0.50 

gals 


M200 




ALL 


0.80 ± 0.05 


1.54 ± 0.09 


0.64 ± 0.06 


183 


177 


N 0.50 
gals 


M200 


P 


Sub- free 


0.79 ± 0.06 


1.66 ± 0.12 


0.59 ± 0.07 


119 


117 


T 0.50 
^opt 


M200 




ALL 


0.69 ± 0.07 


1.74 ± 0.13 


0.79 ± 0.08 


183 


181 


t 0.50 
^opt 


M200 


P 


Sub- free 


0.70 ± 0.07 


1.72 ± 0.13 


0.70 ± 0.08 


119 


116 
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Table 7. Fit parameters for the relations connecting M200 and ap with N 9a ; s , L op t and 
Lx ■ These last parameters are determined within R.200- Masses are in units of 10 14 Mq and 
velocity dispersion in km s _1 . The parameters of each relation are listed in column 1. In the 
second column, we indicate the method used for estimating mass or ap. Mass can be derived 
from the M200-TX relation of POP05 or from the caustic (RD06). Velocity dispersion values 
come from the latter (RD06). Relations derived with the first method (masses estimated 
from Tj) are based in the 21 clusters with temperature available. Results obtained with the 
second method (considering the caustic values) use the 53 CIRS clusters. The pivot points 
are the same as listed in the previous tables when considering the CIRS sample. For the 
21 rich clusters with Tj available the pivot points are assumed the same. The remaining 
columns have the same meaning as in the tables above. 



Relation 
X 


Y 


Mass 


A 


B 


cr lnY\X 


Nt ot 




gals 


M200 


X-ray 


0.90 ± 0.10 


0.92 ± 0.10 


0.34 ± 0.08 


21 


20 


T -R200 
^opt 


M200 


X-ray 


0.99 ± 0.10 


0.85 ± 0.11 


0.33 ± 0.08 


21 


20 


L x 200 a 


M200 


X-ray 


1.15 ± 0.14 


0.39 ± 0.10 


0.36 ± 0.11 


20 


20 


Lf 200 6 


M200 


X-ray 


1.02 ± 0.12 


0.47 ± 0.08 


0.25 ± 0.07 


19 


18 


N-R200 
gals 


ap 


Caustic 


6.34 ± 0.03 


0.27 ± 0.06 


0.18 ± 0.02 


53 


52 


t i?200 
^opt 


a P 


Caustic 


6.36 ± 0.03 


0.25 ± 0.05 


0.18 ± 0.03 


53 


52 


L« 2 °°a 


ap 


Caustic 


6.36 ± 0.03 


0.12 ± 0.03 


0.21 ± 0.03 


52 


51 


L« 200 b 


ap 


Caustic 


6.35 ± 0.03 


0.12 ± 0.03 


0.21 ± 0.03 


49 


48 


N-R200 
gals 


M200 


Caustic 


0.64 ± 0.14 


1.35 ± 0.34 


0.77 ± 0.22 


53 


53 


t -R200 
^opt 


M200 


Caustic 


0.72 ± 0.12 


1.26 ± 0.29 


0.74 ± 0.18 


53 


53 


L« 2 °°a 


M200 


Caustic 


0.73 ± 0.13 


0.56 ± 0.19 


0.87 ± 0.12 


52 


52 


L« 200 6 


M200 


Caustic 


0.68 ± 0.12 


0.54 ± 0.15 


0.87 ± 0.12 


49 


49 



a— L^ 200 estimated with the "annulus" background; 6— L^ 200 estimated with the "frame" 
background 
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Table 8. Slope and orthogonal scatter of different scaling rela- 
tions, from this work or from other authors. Masses are in units 
of 10 14 Mq and velocity dispersion in km s — 1 . The parameters 
involved in each relation are listed in column 1; the sample con- 
sidered is in column 2; the source in columns 3; while the slope 
and scatter are in columns 4 and 5, respectively. 



Relation 
X 


Y 


Sample 


Source 


B 


a lnY \X 


MH500 
gals 


M 50 o 


ALL 


This work 


1.04 


0.43 


N-R200 

gals 


M 200 


ALL 


This work 


1.07 


0.46 


t H500 
^opt 


M500 


ALL 


This work 


1.06 


0.44 


t R'2{)0 
^opt 


M 200 


ALL 


This work 


1.11 


0.49 


t Rr>oo 

^opt 


M500 


Extended 


Pop05 


1.09 


0.50 


T i?200 
^opt 


M 200 


Extended 


Pop05 


1.06 


0.51 


L« 5 °°a 


M500 


ALL 


This work 


0.62 


0.47 


L^ 200 a 


M200 


ALL 


This work 


0.67 


0.56 


L-R500J 


M500 


ALL 


This work 


0.59 


0.46 


T R200l 


M200 


ALL 


This work 


0.61 


0.55 


Lx 


M 500 


Extended 


Pop05 


0.59 


0.59 


Lx 


M200 


Extended 


Pop05 


0.60 


0.54 


Lx 


M 




Sta06 


0.63 


0.43 


Lx 


M500 




Vik08 


0.62 


0.25 


Lx 


M 




Ett04 


0.53 


0.25 


t R600 
^opt 


Up 


ALL 


This work 


0.55 


0.24 


T H200 
^opt 


dp 


ALL 


This work 


0.56 


0.21 


T H500 
^opt 




Extended 


Pop05 


0.43 


0.17 


t K200 
^opt 


Up 


Extended 


Pop05 


0.43 


0.17 


L« 500 a 


ap 


ALL 


This work 


0.26 


0.29 


L H500j, 




ALL 


This work 


0.24 


0.28 


Lx 


ap 


Extended 


Pop05 


0.27 


0.25 


a- L« 20(1 


estimated with the " 


annulus" background; 


6- L« 200 



estimated with the "frame" background. The w or ks of Pop05, 
Sta06 , Vik08 and Ett04 are iPopesso etaLl ||2005|1; IStanek et alj 
[|2006h : IVikhlinin et alj [|2008h and lEttori etaD d2004h re- 
spectively. Note that the results from IStanek et alj ||2006|) are 
dependent of the cosmology. When considering results from the 
third year of WM AP, they find a red uced scatter c;„ A f |L X = 0.25, 
consistent to IVikhlinin et al 
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Table 9. Fit parameters for the relations between X-ray luminosity (Lx) and optical properties (N ga i s 
and L op t) within three different apertures. The aperture used for computing Lx is always the same as 
indicated for N ga ; s or L op t. Lx is measured in units of 10 44 erg/s. The background type considered for 
the Lx estimates is listed in the second column. The sample used for deriving the relations is shown 
in the third column. The remaining columns are as in Table 1. The pivot points are C = 25 when 



considering 


; richness, and C = 


0.40 10 12 L Q 


for L op i. 












Relation 
X 


Y 


Dl,„ 1"% 

LiKg lype 


Sample 


A 


B 




a lnY\X N «ot 






N 0.50 
gals 


Lx 


Annulus 


ALL 


-1.51 ± 0.09 


2.28 ± 


0.20 


0.83 ± 0.13 


105 


105 


N-R500 
gals 


Lx 


Annulus 


ALL 


-1.69 ± 0.09 


1.64 ± 


0.10 


0.62 ± 0.09 


104 


103 


N-R200 
gals 


Lx 


Annulus 


ALL 


-1.85 ± 0.12 


1.59 ± 


0.12 


0.63 ± 0.09 


97 


97 


T 0.50 
U opt 


Lx 


Annulus 


ALL 


-1.68 ± 0.12 


2.38 ± 


0.23 


0.97 ± 0.19 


105 


104 


T -R500 
U opt 


Lx 


Annulus 


ALL 


-1.77 ± 0.10 


1.62 ± 


0.11 


0.65 ± 0.10 


104 


103 


t -R200 
U opt 


Lx 


Annulus 


ALL 


-1.84 ± 0.13 


1.49 ± 


0.12 


0.64 ± 0.10 


97 


97 


N 0.50 
gals 


Lx 


Frame 


ALL 


-1.47 ± 0.09 


2.24 ± 


0.20 


0.82 ± 0.13 


103 


103 


N-R500 
gals 


Lx 


Frame 


ALL 


-1.71 ± 0.09 


1.69 ± 


0.11 


0.68 ± 0.10 


98 


98 


N-R200 
gals 


Lx 


Frame 


ALL 


-1.88 ± 0.12 


1.62 ± 


0.12 


0.75 ± 0.11 


90 


90 


T 0.50 


Lx 


Frame 


ALL 


-1.63 ± 0.12 


2.34 ± 


0.23 


0.96 ± 0.20 


103 


102 


T -R500 
U opt 


Lx 


Frame 


ALL 


-1.73 ± 0.10 


1.61 ± 


0.11 


0.68 ± 0.11 


98 


97 


t -R200 
^opt 


Lx 


Frame 


ALL 


-1.99 ± 0.13 


1.59 ± 


0.13 


0.78 ± 0.14 


90 


89 



Table 10. Fit parameters for the relations between X-ray temperature Tx (in keV) and 
optical properties (N ga i s and L op t) within three different apertures. Only the 21 clusters 
with x-ray temperature available in BAX are used in the fits. That is indicated in the 
second column ("sample"). The remaining columns are as in the previous table. The pivot 
points arc C = 60 when considering richness, and C = 1.10 10 12 Lq for L op t- 



Relation 
X 


Y 


Sample 


A 


B 


< J lnY\X 


N tot 




N 0.50 

gats 


Tx 


BAX 


1.52 ± 0.05 


0.73 ± 0.11 


0.20 ± 0.05 


21 


21 


N-R500 
gals 


Tx 


BAX 


1.22 ± 0.06 


0.50 ± 0.07 


0.19 ± 0.04 


21 


21 


N-R200 
gals 


Tx 


BAX 


1.05 ± 0.06 


0.56 ± 0.06 


0.21 ± 0.05 


21 


20 


T 0.50 
L opt 


Tx 


BAX 


1.56 ± 0.07 


0.69 ± 0.13 


0.23 ± 0.05 


21 


21 


T -R500 
U opt 


Tx 


BAX 


1.21 ± 0.06 


0.52 ± 0.07 


0.20 ± 0.04 


21 


20 


t .R200 
u opt 


Tx 


BAX 


1.10 ± 0.06 


0.51 ± 0.06 


0.21 ± 0.04 


21 


20 



NoSOCS in SDSS. II 



Table 11. Fit parameters for the Rsoo-N 3a ; s , R.200-N 9a ; s , Rsoo-Lopt and R200-L O pt relations 
using three different apertures for computing N 9a ; s or L op t- The physical radii are measured 
in Mpc. The remaining columns are as in the previous table. The pivot points arc C = 25 
when considering richness, and C = 0.40 10 12 Lq for L op t- 



Relation 
X 


Y 


Sample 


A 


B 


G lnY\X 


N tot 




N 0.50 
gals 


R500 


ALL 


-0.12 ± 0.02 


0.47 ± 0.02 


0.19 ± 0.01 


180 


173 


N-R500 
gals 


R500 


ALL 


-0.19 ± 0.01 


0.37 ± 0.01 


0.14 ± 0.01 


180 


175 


N 0.50 
gals 


R200 


ALL 


0.20 ± 0.02 


0.48 ± 0.02 


0.19 ± 0.01 


180 


173 


N-R200 
gals 


R200 


ALL 


0.05 ± 0.02 


0.39 ± 0.02 


0.17 ± 0.01 


180 


177 


T 0.50 
^opt 


R500 


ALL 


-0.16 ± 0.02 


0.56 ± 0.04 


0.24 ± 0.02 


180 


178 


t A500 
U opt 


R500 


ALL 


-0.22 ± 0.01 


0.40 ± 0.02 


0.17 ± 0.01 


180 


175 


t 0.50 
^opt 


R200 


ALL 


0.16 ± 0.02 


0.56 ± 0.04 


0.24 ± 0.02 


180 


178 


t i?200 
^opt 


R200 


ALL 


0.01 ± 0.02 


0.40 ± 0.02 


0.18 ± 0.01 


180 


176 



Table Al. X-ray luminosity and temperature of the 183 NoSOCS plus CIRS clusters. 



name 


Lx (annulus) 


Lx (frame) 


L x (box) 


T x (BAX) 


T x (interp. RASS) 


NOTE 




(10 44 crg s" 1 ) 


(10 44 crg s" 1 ) 


(10 44 crg s" 1 ) 


(keV) 


(keV) 




NSCS J121847+484410 


0.049 ± 0.019 


0.052 ± 0.020 


0.051 ± 0.019 




0.7 


SD 


NSCS J011502+002441 


0.404 ± 0.052 


0.497 ± 0.058 


0.497 ± 0.056 


r, „ +0.09 

-0.09 


2.1 


SD 


NSCS J100242+324218 


0.395 ± 0.048 


0.342 ± 0.045 


0.371 ± 0.046 




2.1 


SD 



Note. - A portion of this table is shown here for guidance regarding its form and content. A full version is available in the 
electronic edition of the MNRAS. 



